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ABSTRACT 

We describe a method for computing the biases that systematic signals introduce in 
parameter estimation using a simple extension of the Fisher matrix formalism. This 
allows us to calculate the offset of the best fit parameters relative to the fiducial model, 
in addition to the usual statistical error ellipse. As an application, we study the impact 
that residual systematics in tomographic weak lensing measurements. In particular we 
explore three different types of shape measurement systematics: (i) additive systematic 
with no rcdshift evolution; (ii) additive systematic with redshift evolution; and (iii) 
multiplicative systematic. In each case, we consider a wide range of scale dependence 
and redshift evolution of the systematics signal. For a future DUNE-like full sky sur- 
vey, we find that, for cases with mild redshift evolution, the variance of the additive 
systematic signal should be kept below 10~ 7 to ensure biases on cosmological parame- 
ters that are sub-dominant to the statistical errors. For the multiplicative systematics, 
which depends on the lensing signal, we find the multiplicative calibration mo needs 
to be controlled to an accuracy better than 10 -3 . We find that the impact of system- 
atics can be underestimated if their assumes redshift dependence is too simplistic. We 
provide simple scaling relations to extend these requirements to any survey geometry 
and discuss the impact of our results for current and future weak lensing surveys. 
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1 INTRODUCTION 

Weak gravitational lensing, or 'cosmic shear' is undergoing 
a phase of rapid expa nsion (see lRefregie r 2003; M unshi et al.l 
2006: lHoekstrall2003l . for reviews) with many future surveys 
and instruments being planned (e.g. DUNIq PanSTARRS3, 
DEifl SNAlQ and LSST0). Central to the planning and de- 
signing of these instruments is our ability to predict the 
uncertainties that such measurements will achieve on the 
cosmological parameters. To this end, the Fisher matrix has 
become a widely used tool in cosmology for calculating their 
covariance matrix. However, a limitation of this approach is 
that it is only able to account for statistical errors, i.e. ones 
that cause an enlargement of the error bars, and is not well- 
suited for treatment of systematic errors, which can intro- 
duce biases that move the measured central value relative 
to its true value. One approach that is commonly taken to 
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overcome this limitation is to treat the systematic errors in 
the same way as statistical errors and to marginalise over 
possible values. This introduction of nuisance parameters, 
in addition to the cosmological parameters, causes the er- 
ror ellipses to expand. A more accurate approach, which we 
use here, is to directly calculate the bias that the system- 
atic signals will introduce. This bias will tend to offset the 
central value to the measurements from the true values, as 
shown in figure [1] Since the computations needed for this 
calculation are very similar to those performed in the stan- 
dard Fisher matrix analysis, extending the current Fisher 
matrix analysis to include a calculation of bias is relatively 
straightforward. 

We apply this formalism to study the impact of resid- 
ual systematics on tomographic cosmic shear surveys. In 
particular, we consider systematics arising in the measure- 
ment of galaxy shapes after correction of instrumental ef- 
fects (such as the Point Spread Function). We consider both 
additive and multiplicative systematics and explore a wide 
range of scale and re dshift dependences. In an earlier work, 
iHuterer et al.l ([2006) considered the impact of photometric 
calibration errors and power law shape systematics using the 
Fisher matrix formalism. A similar bias formalism was intro- 
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duced by iHuterer and Takadal i|2005l ) and applied to theo- 
retical uncertainties in modeling the matter power spectrum 
with N-body simulations. Our work expands upon these ear- 
lier works, by appying the bias formalism to a broad set of 
shape systematics and by studying the joint impact of sys- 
tematic and statistical errors in current and future surveys. 

This paper is organised as follows. In section 2, we de- 
scribe the formalism that we use to quantify systematic bi- 
ases. In section 3, we apply our formalism to cosmic shear 
surveys by exploring the effect of three types of shape mea- 
surement systematics: (i) additive with no redshift evolu- 
tion; (ii) additive with redshift evolution; and (iii) multi- 
plicative. For each type, we consider several possibilities for 
their scale dependence: (i) log-linear systematics; (ii) sys- 
tematics that have the same shape as the lensing signal; and 
(iii) systematics that mimic a small change in the cosmolog- 
ical parameters. In section 4, we study the impact of the 
systematics in the design of future surveys. Our conclusions 
are summarised in section 5. 



2 GENERAL BIAS FORMALISM 

In a weak lensing survey, the observed power spectrum de- 
rived from the shapes of background galaxies is given by 



Clens . y—isys . ^noisc 



(i) 



where each term corresponds to, respectively, the lensing 
signal (Cj cns ), residual systematics (C| ys ), and noise arising 
from measurement errors and intrinsic shape noise (C7° lse ). 
An estimator of the weak lensing shear power spectrum can 
thus be defined as 
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where it is assumed that the residual systematics is unknown 
and therefore uncorrected. The errors of this estimator are 
given by 
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where / s k y is the fraction of the sky covered by the survey. 

The measurement of this power spectrum can then be 
used to constrain a set of cosmological parameters pt. For 
this, we form the usual statistic 
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An estimator for the parameters pi is then defined such that 
dX 2 (Pi)/dpi = 0. 

Neglecting the dependence of the errors ACe in the pa- 
rameters, the covariance matrix of the parameters 

cov]pi,pj] = {{pi - (pi)){pj - (£,})) = (F _1 )ij (5) 

is then given by the inverse of the Fisher matrix 

F a = A<7 i" 2 ~t~ ~tr ■ («) 



dpi dpj 

It is also easy to show that, for small residual systematics, 
the bias of the parameter estimator is given by 



b[Pi] = {Pi) - ( P T UC ) = (F- 1 ) ij B j 



where p* rue is the true value of the parameters, the summa- 
tion conv 
given by 



B 



, dCf 
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This simple expression is similar to that for the Fisher 
matrix and is therefore a convenient way to evaluate 
the impact of residual sy stem atics and has also u sed by 
IHuterer and Takadal 1120051 ) and IHuterer et all l|2006ft . Note 
that this expression, while derived for the measurement of 
a single weak lensing power spectrum, is general and can be 
applied to any estimation of a model using a \ 2 fit m the 
presence of residual (unknown) systematics. In particular, it 
can be easily generalised for the case where several power 
spectra are considered, such as in weak lensing tomography, 
and if the power spectrum estimators at different multipoles 
are correlated. 

Figure [T] illustrates the principles of our formalism for 
the case of the measurement of dark energy equation of 
state parameters (see detailed discussion below). The black 
dashed ellipse shows the statistical error estimates for these 
parameters derived using the fisher matrix. This error el- 
lipse is centred on the fiducial model. The red error ellipse 
includes the systematics bias discussed in equation [8] We 
see that systematics have an additional effect of shifting the 
centre of the error ellipses away from the fiducial model. 

In general, it is convenient to consider the total error 
covariance matrix which considers deviations from the true 
value of the parameters and is given by 

tcov[pi,&] = ((pi-pT ne )(pj- pf™)) = (•F _1 )y+&[p<]&[Pj]-(9) 

This includes both statistical and systematic errors, as op- 
posed to the the standard covariance matrix (Eq. [5} which 
only considers deviations from the mean values of param- 
eter estimators and thus only includes statistical errors. In 
particular, the diagonal elements of the total error matrix 
give the Mean Square Error (MSE) of a parameter pt 



MSE[p l ]=a 2 [p 1 ]+b 2 [p i 



(10) 



where the statistical error variance is a 2 [j5i] = (F~ 1 )a. It is 
this total error which needs to be minimised when optimis- 
ing future surveys rather than the statistical error alone. An 
interesting criterion is to define a tolerance on the system- 
atics such that they do not dominate over statistical error. 
This is verified when 



b\pi] < a\ 



(11) 



(7) 



for all or a selected subset of the paramters pi. In the fol- 
lowing, we will apply this formalism to cosmic shear and 
derived the systematics tolerance for future surveys. 



3 SHAPE SYSTEMATICS FOR WEAK 
LENSING TOMOGRAPHY 

Weak lensing tomography places constraints on cosmologi- 
cal parameters by measuring the statistics of the shear field 
as a function of redshift. This is done in practice by di- 
viding galaxies into redshift bins and measuring both the 
auto-correlation of the shear signal within bins as well as 
the cross-correlation between bins. The shear signal itself is 
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Figure 1. Illustration of the distinction between statistical errors 
which can be estimated using the Fisher matrix and the biasing 
effect of residual systematics, which can be estimated using the 
present formalism. The black dashed line shows the results of 
a Fisher matrix calculation with the cross showing the fiducial 
model that has been used. The solid red line shows the error 
ellipses when the effects of a systematic signal are also included. 
We see that the systematic errors can induce a bias that moves 
the central value relative to the fiducial model. 



measured using the shapes of distant galaxies. For this pur- 
pose, the observed galaxy images must be corrected from 
the effects of the point spread function (PSF) which is mon- 
itored using stars in the image. Errors in the PSF decon- 
volution and galaxy shape estimation induces errors in the 
measured galaxy shapes, and hence lead to errors in shears. 
These shape measurement errors generally have spatial cor- 
relations (since the PSF itself has a spatial variation) which 
can be described by its power spectrum C| ys . In this work, 
we decompose residual systematic power spectrum into two 
parts: an additive term (i.e. independent of the lensing sig- 
nal) and a multiplicative term (i.e. dependent on the lensing 
signal). Each term and its impact on the systematics power 
spectrum is described in AppendbtAl 

In the following, we consider a flat cosmological model 
with 7 parameters listed in table [1] In particular, the 
evolution of the dark energy equation of state parameter 
is assumed to take the form w(a) = w n + (a n — a)w a , 
with wo, listed in the table, corresponding to a pivot 
point of a n = 1. For the cosmic shear survey parameters, 
we focus on the DUNE-like 'shallow' survey described in 
lAmara and Refregierl |2006), namely a 20,000 sq. degree 
survey containing 35 galaxies per amin 2 with ct 7 = 0.25 
and Zm = 0.9, with the galaxies divided into 5 redshift 
bins so that each bin contains the same number of galax- 
ies .^We_jssmne that t he overall galaxy distribution is given 
bv lSmail et ail (| 19941 ). 



P{z) 



z exp 



where we set a = 2 and = 1.5. The median redshift of the 
survey, z m , is then used to set zq ~ z m /1.412. We use multi- 
poles in the range 10 < £ < 2 x 10 4 . The lensing tomo g raphy 
for malism we use is also descri bed in iHu and Jainl (|2004h 
and lAmara and Refregierl (120061 ). Table [1] shows the central 
values of our fiducial model and the marginalised Fisher ma- 
trix errors on each of the cosmological parameters for this 
survey. These errors come from the same calculation used to 
produce the results (dashed curve) shown in figure [T] The 
errors on the dark energy parameters are sometimes stated 
in terms of w n corresponding to the pivot point a n where 
w n and w a are uncorrelated. For our survey the marginalised 
error on w n is 0.02. 



3.1 Additive Term 

First, we consider an additive systematic signal, which, by 
definition is not correlated with the lensing signal. Such a 
systematic could, for instance, result from residual errors in 
the correction of the PSF. This signal can in general have 
an arbitrary scale dependence and may also depend on the 
galaxy redshift. 

To quantify the amplitude of the systematic signal it is 
convenient to consider its variance 



1 

2^ 



\Cf s \£(£ + l)dln£, 



(13) 



where the absolute value sign is included to account for pos- 
sible changes of sign. To explore the possible scale depen- 
dence of the systematic signal, we consider the following 
three classes of shape of the systematic power spectrum: 

• The first is a log-linear systematic: 

£(£ + l)Cf s = Ao(n\og w {£/£ ) + l), (14) 

where £o is a reference scale, n a scaling parameter and 
Ao is a normalisation. This parametrisation allows for the 
possibility that the residual power spectrum of systematic 
signal, after correction, be positive or negative and may 
transition from one to the other. More specifically, the 
quantity £(£ + l)C| ys scales linearly with log^ and goes 
through the point £(£ + l)C| ys = Ao at £ — £o and has 
a slope of nAo.In general, each correlation function could 
have its own normalisation (i.e. Aij, rather than Ao). For 
instance this would be the case if the systematic signal had 
a redshift dependence. 

• Secondly, we explore systematic signals that have the 
same shape as one of the lensing power spectra, C\ cnB . Since, 
in this work, the galaxies have been divided into 5 redshift 
bins, this gives us 15 possible power spectra to investigate: 



Cf = Aid 



(15) 



(12) 



where Cj J is the correlation spectrum between bins i and 
j, and Ai is a normalisation. Once again each correlation 
function could have a different normalisation. For instance 
this would be the case if the systematic signal had a redshift 
dependence. 

• Finally, we explore the systematic shape that should 
have the greatest impact on our measured cosmological pa- 
rameter estimation, namely systematics that exactly mimic 
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Parameter Q m wq w a h erg Qf, n 

Central value 0.28 -0.95 0.00 0.72 1.0 0.046 1.0 

Marginalised errors 0.006 0.06 0.19 17 0.008 4 0.009 

Table 1. Cosmology parameters for our fiducial model with marginalised Fisher matrix errors for the survey we consider. 
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Figure 2. Different classes of additive systematic signals we con- 
sider. In each case, the systematics power spectra shown have 
been normalised so that they induce a a bias on wq equal to 
the systematics error (b[uio] = <r[jiio] = 0.06). The top left panel 
shows examples of the log-linear models for n=-1.4, -1.0, -0.6, - 
0.2, 0.2 and 0.6; the top right panel shows systematics that have 
the same shape as the lensing signal (all 15 powerspectra shapes 
are shown); and the bottom left panel shows signal that have the 
shape of dCf 3 / dp a , where C| 3 is the auto-correlation power spec- 
trum of the third redshift slice and p a are the parameters in our 
fiducial model. Finally, the bottom right panel shows systematic 
signals that have the same functional form as dCl 1 /dwo , where 
CV is the lensing correlation function between bins i and j . 



the effect of a change in one of the cosmological parameters 
(Pa): 



C7 S = A: 



(16) 



whereA2 is another normalisation factor. 



Figure [2] shows examples of these different classes of 
systematic signals. Each of the systematic signals have been 
normalised so that they introduce a bias on the equation of 
state parameter wo of 0.06, which is the same level as the 
marginalised statistical error on this parameter calculated 
from a standard Fisher matrix approach see Table [T] 



3.1.1 Additive term with no redshift evolution 

We first consider additive systematics whose amplitude is 
independent of redshift. Figure [3] shows the bias introduced 
to each parameter in our log- linear model (Eq. [14] and top 



Figure 3. The bias of all 7 parameters considered in our fiducial 
model produced by the log-linear systematic as a function of n, 
the scaling in equation 1141 Aq is chosen so that ctJ = 4 X 10~ 7 . 



left panel of figure [2j|) as function of the scaling parameter 
n. The 7 curves in the figure show the bias for a system- 
atic signal with an amplitude such that CTg ys = 4 x 10 -7 . As 
can be seen from equations 151 and IT3l the bias scales linearly 
with 0"g ys . The level of the bias depends on the cosmological 
parameter considered and the absolute value of the bias re- 
mains roughly flat with n. The apparent dips around n = — 1 
correspond to a change of sign of the biases (e.g. goes from 
positive to negative) . The black-solid and blue-dashed curves 
show that a log- linear systematic with a^ ys = 4 x 10 -7 will 
introduce a bias of ~ 0.05 on tuo and ~ 0.1 on w a , for most 
values of n. The range of possible bias values is large depend- 
ing on which cosmology parameter is under scrutiny (note 
that the j/-axis of figure [3] is on a log scale). For example, 
the absolute bias induced in f2 m is typically ~ 0.005 for a 
wide range of n values. 

Although the value of the bias on the Q m is smaller than 
that for too, a more relevant quantity is the size of the bias 
relative to the statistical errors (which are shown in table 
[TJ • Figure U shows the ratio of the bias to the marginalised 
Fisher matrix errors of all 7 cosmology parameters for the 
log-linear model with a^ ys — 4 x 10 -7 . Expressed in this 
manner, we see that many of the cosmological parameters 
respond to the same extent to the presence of a systematic 
signal (note that the j/-axis is now on a linear scale) . We also 
see more clearly that as we vary n, the bias on any given 
parameter can change sign. 

The grey shaded region in the figure is where the bias 
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dominates over statistical errors, i.e. where \b(pi)\/a(pi) > 1. 
With this in mind, we can now adjust the level of the the 
systematic (<Ts ya ) and find our tolerance (o" ays _ tol ), which, 
from equation [11] is the maximum value that cr ays can have 
before the potential bias on one of the cosmology param- 
eters breaches our acceptable threshold (i.e. before one of 
the curves in figure [4] enters the grey area). Figure 4 shows 
that, for the log-linear systematic, setting a threshold of 
0"s ys _toi — 4 x 10~ 7 would be sufficient to meet this criterial 
for the dark energy parameters. 

Appendix [B] shows the results for the other classes 
of systematics. Table IB1I shows the ratio of the bias to 
marginalised statistical errors for the residual systematic 
signals shown in the top right panel of figure [2] table IB2I 
correspond to the bottom left panel figure [2] and IB3I are 
the results for the bottom right panel. We see that, if the 
residual systematic has the same shape as the lensing sig- 
nal, then a tolerance level of <r aya _ tol = 4 x 1CT 7 is suffi- 
cient to keep \b(pi)/a(pi)\ < 1 for the wo and w a param- 
eters. However, table IB3I shows that a systematic level of 
o"g ys _ tol = 4 x 10 -7 can lead to a ratio greater than one for 
wo for some systematic power spectrum shapes, for instance 
Cfys oc dCf3/dwo- It is also interesting to note the the most 
stringent constraints on systematics does not always come 
from dark energy. For instance for a systematic that has a 
shape C^ ys oc dCf 5 /dwo, it is the constraint from Q, m that 
would dominate. 

Focusing again on dark energy parameters, figure [5] 
shows the biases (as a fraction of statistical errors) for wo 
and w a for all the additive systematics we have considered. 
Once again the grey area shows regions where the bias would 
dominate over the statistical errors. We see that a system- 
atic with o"s ys _ tol = 4 x 10 -7 is not always sufficient (this 
can also be seen in tables [B2l and IB3T1 . Instead a more com- 
fortable tolerance level is a^ ys _ tol = 3 x 1CP 7 , which we can 
be considered as the target for the control of additive sys- 
tematics with no redshift evolution. 

3.1.2 Additive term with redshift evolution 

The additive term discussed in the previous section was for a 
systematic signal that is the same for all tomographic power 
spectra, i.e. independent of the galaxy redshifts. An exten- 
sion of the above results is to consider an additive term that 
evolves with redshift. This, for example, could be due to 
the fact that more distant galaxies tend to be smaller than 
nearby galaxies, making their shape measurements more 
prone to errors in PSF deconvolution. We can model such an 
effect by introducing a simple redshift scaling to the additive 
part of equation I A5I 

add add r -t , \0 n f -\ n\ 

7 =7o (l + 2m) Pa , (17) 

where z m is the median redshift of a tomographic redshift 
bin. The systematic power spectrum now depends on the 
redshift of the two galaxies being correlated and, for sim- 
plicity, we assume that this redshift dependence can be sep- 
arated from the angular dependence such that, 

c^izLzl) = (i+4 l )""(n-a4)""Cf(4 l = o,*£ = o), (is) 

where z % m and z 3 m are the median redshifts of bins i and j, 
and Cf s (z % = 0, z 3 — 0) is the systematic auto-correlation 
power spectrum of galaxies at redshift z=0. 




-2 f- 1 

-3 -2-10 1 

n 



Figure 4. The ratio of the systematic bias to the marginalised 
statistical errors for each of the cosmological parameters as a func- 
tion of n, the scaling parameter of the log-linear scaling model. 
The amplitude has been set to o- aya = 4 x 10 — 7 . The grey area 
corresponds to regions where the bias dominates over the statis- 
tical errors (i.e. \b(pi)/cr(pi)\ > 1). We see that it is not always 
the dark energy equation of state parameters that drive the sys- 
tematic requirements. 
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Figure 5. Impact of additive systematics (with no redshift evo- 
lution) on the dark energy parameters, u>o and w a , for all classes 
of scale dependence considered. As before, the grey area corre- 
sponds regions where the bias dominates over the (marginalised) 
statistical errors. In all cases, we set CT ays = 4 x 10 — 7 . 
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Figure 6. The ratio of bias to statistical errors for additive sys- 
tcmatics with redshift evolution for the different cosmological pa- 
rameters. The shape of the systematics has been chosen to be 
Cg ys oc dCf 3 /cLwq and j3 a is the power-law exponent describing 
the redshift evolution (see text). In all cases, the systematic auto- 
correlation variance at z = is set to cr^ yB = 10 — 7 . Note that this 
is a systematic amplitude that is a factor of 4 smaller than that 
used in earlier figures. 



Figure [6] shows the bias (as a fraction of marginalised 
errors) for different values of /3 a and for systematic shape of 
the form C s t vs oc dCf 3 /dwo (one of our worst-case systematic 
shapes) and an amplitude such that (7s ys = 1 x 10~ 7 . Once 
again we see that it is not always the dark energy equation 
of state parameters that drive the systematic requirements, 
and that Sl m and as can be as restrictive as the two w pa- 
rameters. We also see large values of f3 a cause the bias to 
become very large relative to the statistical errors. For the 
dark energy parameters, we therefore place a requirement of 



I3 a < 1.5 and cr tol (z 



0. 



o) < i x lcr 



3.2 Multiplicative Term 



Having investigated the effect of additive systematics, we 
now investigate the case of multiplicative systematics which 
are correlated to the lensing signal. This is most easily done 
by studying the impact of different values of the multiplica- 
tive calibration m defined in equation IA5I and IA8I First, 
we consider a simple scaling of m with the median redshift, 
since the redshift can have an impact on galaxy properties, 
such as size and magnitude, which in turn have an impact 
on the accuracy with which their shapes are measured. The 
redshift scaling is thus taken to be 



mo(l + z m ) 



(19) 



Figure [7J shows the fractional bias as a function of f3 m 
with mo = 5 x 10 -3 . If there is no redshift dependence (i.e. 
if p m = 0) then setting a tolerance of mo < m to ; = 5 x 1CP 3 
is more than sufficient to keep the biases at a subdominant 



levels. Once again, it is not the equation of state parame- 
ters that place the most stringent constraints on systemat- 
ics. In particular, the primordial shape parameter n is the 
most sensitive. If one is only concerned with dark energy 
parameters, we can set a limit of /3 m < 1.5 for this value 
of mtoi- Figure [7] also shows some distinct features that are 
not apparent in figure [6] namely that the biases change sign 
as /3 m varies. Since /3 m controls the relative amplitude of 
the systematics as a function of redshift, this suggests that 
not all the power spectra are affected in the same way. For 
instance, for wo, when the error on the low redshift cor- 
relation functions dominates (i.e. for small /3 m ) the bias is 
negative, whereas when the errors on the high redshift cor- 
relation functions dominates (i.e. for large /3 m ) the bias is 
positive. This is important since it shows that there exists 
a trade-off between high and low redshift biases. Moreover, 
this indicates that if m transitions through zero at some 
redshift, the multiplicateive systematics are likely to have a 
greater effect. This is analogous to what we have found for 
the additive term. 

From figure [7] we would thus expect that if we are able 
to constrain mo to 5 x 10 -3 with [3 m — 1, we would ex- 
pect that the bias on the dark energy parameters would be 
significantly subdominant to the statical errors. Figure [5] 
however, tells a different story. For any experiment the tar- 
get is for m to be as close to zero as possible. Since m can 
be either positive or negative (or infact transition from one 
to the other) one can image defining a tolerance envelope 
such that the unknown systematic can be anywhere in this 
region. This is illustrated in figure [5] The shaded area shows 
the uncertainty envelope enclosed by m = imo(a + z)^ m , 
with mo = 5 x 10 -3 and f3 m = 1. In blue we show the simple 
model used from equation 1191 where m is always positive 
and the red dashed curve, which was used to produce figure 
[8] is for a systematic that has the form, 



= ^-atan^a m (z m - z T )^ m (l + z™)' 5 " 



(20) 



where zt is the redshift at which the systematic transitions 
through zero and a m controls the rate of the transition. 
Figure [TU] shows results for a number of a m and /3 m values 
as a function of the transition redshift. From this we are 
able to set the requirements that when we allow m to have 
a zero transition, the tolerance is m -toi = 1 x 10~ 3 for 
Pm < 1.5. 



4 IMPACT ON WEAK LENSING SURVEYS 



In [A mara and Refregieil (2006) we studied how the geome- 
tries and configurations of future lensing surveys can be op- 
timised to derive the strongest cosmological constraints. The 
present results can be used to extend this work, which solely 
included statistical errors, to also include the impact of sys- 
tematic biases. For this purpose, it is useful to consider the 
total mean square error of the cosmological parameters pi, 
which we have defined in equation 1101 

Figure [TT1 shows the root mean square error as a func- 
tion of survey area, with all other properties of the survey 
set to the values listed in section [3] Also shown are the 
impact of different levels of additive systematics with am- 
plitudes of <Tsy S = 10 -5 , 10 -6 and 10 -7 , from top to bottom 
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Figure 7. The ratio of the bias to statistical error for the mul- 
tiplicative systematics for the power-law redshift evolution. The 
multiplicative calibration parameter is set mo = 5 X 10 — 3 . 



respectively. The systematics were chosen to have a shape 
C s t ys oc dCf 3 /dw n without redshift evolution. The figure 
demonstrates that, when designing future surveys, one needs 
to consider both the statistical power of the survey and the 
systematic floor of the survey. For instance, we see that if 
an instrument was to have systematics with a1 ya = 10 -5 
then there is very little point in surveying an area of sky 
bigger than 20 square degrees. Alternatively, we see that if 
we choose to do a 20,000 square degree survey, we will need 
to control the residual systematic to o1 ys ~ 10~ 7 in order 
to reach the statistical potential of the survey. 

Since the requirements depend on the ratio of the bias 
to statistical errors, the requirements for a specific survey 
configurations depends on the statist i cal p ower of that spe- 
cific survey. In [A mara and Refregieil (2006) we investigated 
the statistical power for measuring dark energy of several 
survey configurations. We then gave a simple scaling rela- 
tion for calculation the dark energy Figure of Merit (FoM) 
as a function of the survey properties. We have found that 
these scaling relations hold well and capture most of the 
relevant physics. In the same spirit we now give the scaling 
relation for the systematic requirements. For a survey that 
has an area A s , a galaxy surface density n g (useful for lens- 
ing) and a median redshift z m , the tolerance on the additive 
systematic is, 



Figure 8. The ratio of bias to statistical errors when m is al- 
lowed to transition between positive and negative. Here we have 
set mo = 5 X 10~ 3 and /3 m = 1. With the simple redshift evolu- 
tion used in figure [7] one would conclude that this level would be 
sufficient for the equation of state parameters. Allowing a more 
complicated redshift evolution can however cause the biases to 
dominate. 
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Figure 9. The shaded area (yellow) shows the tolerance envelope, 
with /3 m = 1 and a— 10. Within this zone the residual systematic 
3 can take any form. In blue we show the systematic described in 
(2£fcjuation[T9]and the red dashed curve shows the systematic from 
equation 1201 



and the requirement on the multiplicative systematic is, 
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Figure 10. Same as Results in figure[9]but for different values of 
/3m, and a m . We have set mo = 10 — 3 , i.e. 5 times more stringent 
than those iinfered from a simple scaling that is always positive. 



Figure 11. The root mean square error (RMSE) on the dark 
energy equation of state parameter wo as a function of survey 
area. This error includes both the effects of statistical errors and 
the bias induced by a systematics with a shape, C| HS tx dCf 3 / 'dwo 
with o 2 sys = 10" 5 , 10" 6 and 10" 7 and (3 a = 1. The red line shows 
the uncertainty with no systematic errors. 



In both cases we have found that the evolution with redshift 
must remain weak (/3 < 1.5). 

To illustrate this scal ing relation we c o nside r its impact 
for current day surveys. iBeniamin et al.l ((2007) have per- 
formed a joint analysis of the latest cosmic shear surveys 
together covering about 100 deg 2 with a median redshift 
of 0.78 and roughly 10 galaxies per arcmin 2 on average. 
From our scaling relation we see that, for this analysis, the 
shear calibration m needs to be determined to a precision 
of better than 3%. As we have already stated, this is close 
to the limit of current shear measurement methods in sim- 
ulated, therefore somewhat idealised, conditions. This un- 
derscores the fact that, for future larger surveys, control of 
systematic signals needs to be a priority since we are al- 
ready at the point where they are likely to dominate. With 
the progress currently being made in the STEP program 
it is reasonable to expect the shape measurements meth- 
ods to improve their accuracy for multiplicative systemat- 
ics, and hence be able to keep up with ever more ambitious 
surveys. For the addit ive requirement we find that for the 
IBeniamin et al. (2007) study we need a\ ys < 3 x 10 6 and 
corresponds rou ghly to an additive er r or of c ~ 0.002 (as d e- 
fined in STEP (jrlevmans et all 120061 : iMassev et al.ll2007h ). 
While less attention has been given to additive systemat- 
ics, STEP studies have shown this level is also achievable 
with the best current measurement methods, assuming that 
the measurements are not limited by PSF modeling and in- 
terpolation. A detailed study of requirements for additive 
systematics and PSF calibration will be presented in a later 
paper (Paulin-Henriksson et al in prep.). 



5 CONCLUSIONS 

In this paper, we have outlined a method for computing the 
biases that residual systematics introduce. This approach 
involves a simple extension of the Fisher matrix formalism 
that is now widely used in cosmology to make error forecasts. 
As an application, we have used it to study the impact that 
residual systematic signals will have on future tomographic 
cosmic shear measurements. Specifically, we have explored 
three different types of shape systematic signal affecting to- 
mographic shear power spectra: (i) additive systematics with 
no redshift evolution; (ii) additive systematics with redshift 
evolution; and (iii) multiplicative systematics. The require- 
ment target is then to have all types of systematics close to 
zero. This defines a tolerance envelope for the systematics 
that allows the residual systematics errors in the power spec- 
tra to be positive or negative within its limits. It is impor- 
tant to note that it is the worst systematic possible within 
this limit which drives the requirement, not a marginalised 
average over all systematics. To this end we have investi- 
gated a wide class of possible systematic shapes and used the 
most constraining ones to set our systematic requirements. 
In doing this, we have found that, for both the additive and 
multiplicative parts, it is vital to consider systematics that 
have positive and negative power spectra C| as . For instance 
we see, in the multiplicative case, that investigating only 
power-law behaviour for its redshift evolution (i.e. m is al- 
ways positive) can lead to a factor of 5 underestimation of 
the impact of a systematic within a given tolerance win- 
dow. From our calculation we are able to set the following 
requirements on the survey we have considered (a DUNE- 
like survey covering 20,000 sq. degrees with 35 galaxies per 
arcmin -2 and a median redshift of z m — 0.9): 
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• For both the additive and the multiplicative signals, 
the redshift evolution needs to be weak /3 < 1.5, where the 
errors for a given galaxy scale as (1 + z) 13 . This is not a 
trivial requirement since the shapes of more distant galaxies 
are harder to measure since they are smaller and fainter. 

• The power spectrum of the residual additive shear error, 
that is the part that is not correlated to the lensing signal, 
must be controlled such that its amplitude is al ys < 1CP 7 
(as defined in equation 1 13p 

• The multiplicative part needs to be controlled to a pre- 
cision of mo < 10 -3 , where mo is the shear calibration error. 
This means that we need to be able calibrate shears to an 
accuracy of 0.1%, which is about one order of magnitude 
better than the current best measurement methods are able 
to achieve, as determined by the latest S TEP simulations 
jHevmans et alj|2006l ; IMassev et ai1l2007l) . 

These specific requirements apply to our fiducial survey, 
but we provide scaling relations (Eqs l21l and l22|l which show 
the requirements for any survey geometry. We have shown 
that for current survey covering ~ 100 deg 2 (|Beniamin et al.l 
120071 ). we need a% s < 3 x 10" 6 and m < 0.03. This level 
of accuracy is at the limit of the performance of the cur- 
rent best shear measurement methods, as demonstrated by 
STEP. However, further systematics, such as that arising 
from PSF calibration and interpolation, are not accounted 
for by STEP and can dominate the error budget for future 
surveys. A discussion of the requirements for additive sys- 
tematics and PSF modeling in the context of present and 
future surveys will be presented in a later paper (Paulin- 
Henriksson et al., 2007, in prep). 



APPENDIX A: SHEAR SYSTEMATIC ERROR 
AND SYSTEMATIC POWER SPECTRUM 

The shear at a position can be written as a complex number, 

j(6,z)=- fl (0,z)+ij2(O,z). (Al) 

This shear field can be measured and used to construct a 
set of correlation functions between the shear measured at 
(6,z) and (9',z'), where 9 and 8' are the angular positions 
of two galaxies on the sky with 9' — (9 + <j>), and z, z' are 
their respective redshift. The two correlation functions 
and £_ can then be constructed as, 

t+(4>,z,z') = ( 7l (0,z) 7l (0V)> + { l2 (9,z) l2 (e',z')), (A2) 
and 

Z-(4>,z,z') = ( 7 i(M7i(0',z')> - M0,zh 2 (e', z')). (A3) 

These real space correlation functions can be combined to 
calculate the power spectrum between galaxies at a redshifts 
of z and z' , 

C(£,z,z') = J ^d<f>^(ct>,z,z')J (£cj>)-^(4>,z,z')J A (e<P)),(AA) 

where Jo and J a are the zeroth and the fourth order Bessel 
functions. In weak lensing tomography, where galaxies are 
divided into redshift slices, we use the notation CV , the 
correlation function between slices at redshift z % and z 3 , for 
C{t,z,z'). 

We assume that the observed shear (without noise) is a 
combination of shear due to lensing plus a residual system- 
atic, 7° bs = -y lms + "f sys , an d that the lensing systematics 
can be decomposed into an additive error and a multiplica- 
tive error, 

~1 sys {9, z) = -y add (9, z) + m(z)j len3 (9, z). (A5) 

This is similar to the decompositio n used by the STE P 
collaboration l|Hevmans et al.l 120061 ; IMassev et al.l I2007T I. 
Putting this systematic into equation IA4I we see that the 
systematic power spectrum is composed of two parts, an 
additive part (JJ add ) and a multiplicative powerspectrum 
(C™'), 

c sys (£, z, z) = c add (£, z, z) + C mul {£, z, z'), (A6) 

where the additive power spectrum can be calculated di- 
rectly from the correlation function of ■y add y 

C add (£,z,z') = y n 0#(ef d (^^, 2 ')^oW)-C- dd (0^^V4(^))(A7) 

and the multiplicative part is dependent on the lensing 
power spectrum, 

C mul {£,z,z') = (m(z) +m(z')^C Lms (£,z,z'). (A8) 
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APPENDIX B: TABULATED BIAS RESULTS 

In this section we give the results for the bias calculation of 
several examples of systematics that have the same shape as 
the lensing signal and systematics that resemble dCi/dp. As 
with the results shown in figure [3] we calculate the bias in 
each parameter for systematic signals that have an an am- 
plitude given by (j^ ys = 10~ 7 . Table iBll shows the results for 
systematic in equation IA4I fi.e. the same shape as the power 
spectrum shown in the top right panel of figure [5}, table [B2l 
are for systematics in equation[l6]where we specifically look 
at changes in the power spectrum corresponding to small 
changes in wq. These systematics are shown in the bottom 
right panel of figure [2] Finally, table IB3I shows other exam- 
ples of systematics that are described by 1161 In this case 
these are the systematics shown in the bottom left panel on 
figure [2] 
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Table Bl. The ratio of bias to marginalised statistical errors for each parameter from a systematic that has the same shape as the 
lensing correlation function and cr^ ys = 4 X 10 -7 . These are the systematic signals shown in the top right panel of figure [2] 
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Table B2. The ratio of bias to Fisher matrix errors for each parameter from a systematic that is similar to the change in Cf caused by 
a small change in each of the cosmological parameters. The systematic correlations all have cr^ yB = 4 X 10~ 7 . These are the systematic 
signals shown in the bottom left panel of figure [2] 
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Table B3. The ratio of bias to Fisher matrix errors for each parameter from a systematic that is similar to the change in C?£ corresponding 
to a small change in wo. As before, we use <r^ ys = 4 X 10 — 7 . These are the systematic signals shown in the bottom right panel of figure 
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